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Abstract. We report on the first firm detection of pulsed 7- 
ray emission from PSR B 1509-58 in the 0.75-30 MeV energy 
range in CGRO COMPTEL data collected over more than 6 
years. The modulation significance in the 0.75-30 MeV pulse- 
phase distribution is 5.4<t and the lightcurve is similar to the 
lightcurves found earlier between 0.7 and 700 keV: a single 
broad asymmetric pulse reaching its maximum 0.38 ± 0.03 in 
phase after the radio peak, compared to the offset of 0.30 found 
in the CGRO BATSE soft gamma-ray data, and 0.27 ± 0.01 for 
RXTE (2-16 keV), compatible with ASCA (0.7-2.2 keV). 

Analysis in narrower energy windows shows that the single 
broad pulse is significantly detected up to ~ 10 MeV. Above 
10 MeV we do detect marginally significant (2.1a) modulation 
with an indication for the broad pulse. However, imaging anal- 
ysis shows the presence of a strong 5.6a source at the position 
of the pulsar. To investigate this further, we have also analysed 
contemporaneous CGRO EGRET data (> 30 MeV) collected 
over a nearly 4 year period. In the 30- 1 00 MeV energy window, 
adjacent to the COMPTEL 10-30 MeV range, a 4.4ct source 
can be attributed to PSR B 1509-58. Timing analysis in this en- 
ergy window yields an insignificant signal of 1.1a, but with a 
shape somewhat similar to that of the COMPTEL 10-30 MeV 
lightcurve. Combining the two pulse-phase distributions results 
in a suggestive double-peaked pulsed signal above the back- 
ground level estimated in the spatial analyses, with one broad 
peak near phase 0.38 (aligned with the pulse observed at lower 
energies) and a second narrower peak near phase 0.85, which 
is absent for energies below 10 MeV. The modulation signifi- 
cance is, however, only 2.3a and needs confirmation. 

Spectral analysis based on the excess counts in the broad 
pulse of the lightcurve shows that extrapolation of the OSSE 
power-law spectral fit with index -1.68 describes our data well 
up to 10 MeV. Above 10 MeV the spectrum breaks abruptly. 
The precise location of the break/bend between 10 and 30 MeV 



depends on the interpretation of the structure in the lightcurve 
measured by COMPTEL and EGRET above 10 MeV. 

Such a break in the spectrum of PSR B 1509-58 has recently 
been interpreted in the framework of polar cap models for the 
explanation of gamma-ray pulsars, as a signature of the exotic 
photon splitting process in the strong magnetic field of PSR 
B1509-58. For that interpretation our new spectrum constrains 
the co-latitude to ~ 2°, close to the "classical" radius of the 
polar cap. In the case of an outer-gap scenario, our spectrum 
requires a dominant synchrotron component. 
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1. Introduction 

PSR B1509-58 was discovered as a 150 ms X-ray pulsar in the 
Einstein HRJ and IPC (0.2-4 keV) data from observations per- 
formed inj^79_andj^80_of™pernova remnant (SNR) MSH 
15-52 (Seward & Harnden 1982). The pulsations and the large 
period derivative indicated in the X-ray data were soon con- 
firmed at radio-wavelengths ( Manchester et al. 1982| ), while 
the derived dispersion measure supports its association with 
the SNR. The inferred characteristic age is 1570 year and the 
component of the surface magnetic field perpendicular to the 
spin axis at the magnetic pole is 3.1 x 10 13 GauB, one of the 
highest among the steadily growing sample of radio-pulsars. 
Radio-data collected during an 1 1 yr time span showed that the 
pulsar did not glitch and made a detailed study of its slow-down 



possible (Kaspi et al. 1994). The measured braking index was 
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n = 2.837(1), close to n = 3 expected for a dipole. 

Extensive X-ray studies of PSR B1509-58 and its environ- 
ment have been performed in the early eighties at soft- and 
medium X-ray ener gies using the Einst ein HRI, IPC and SSS 
dSewardet al. 1983j l984) and MPC ( |Weisskopf et al. 1983| ) 
instruments and in the late eighties using the EXOSAT ME 
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and LE instruments (Trussoni et al. 1990). The morphology of 
MSH 15-52 at X-rays is complex with at the north western rim 
of the SNR an excess near the Ha nebula RCW 89 and close 
to the middle of the SNR a clump containing the pulsar sur- 
rounded by a diffuse synchrotron nebula. The spectrum of the 
pulsar is hard with a photon power-law index of ^ — 1.1 ± 0.1 



(Trussoni et al. 199C). The pulse is broad with a duty cycle of 



50% and rather asymmetric: a sharp rise followed by a grad- 
ual decline. At soft X-ray energies (< 2 keV) there is some 
indication that the broad pulse is comp osed of two smaller 
narrowly separ ated pulses (E XOSAT LE, Trussoni et al. 1990 ; 
ROSAT PSPC, |Becker 1994| ). 

More recently, the results from ROSAT PSPC/HRI 
(Greiveldinger et al. 1995, Trussoni et al. 1996, and Bra- 
zier & Becker 



Tamura et al. 1996 



1997) 
and 



and 



ASCA ([Nagase 
Saito et al. 199% 



et al. 



1994, 

observations were 
presented. Using the high-spectral resolution of ASCA in com- 



bination with imaging ( Tamura et al. 1 996 ) and the high-spatial 



resolution of the ROSAT HRI (Brazier & Becker 1997) the 
morphology of the remnant can be explained by the presence 
of several components: the pulsar itself, a non-thermal nebula 
powered by the pulsar with collimated outflow structures (jets) 
and a hot thermal plasma at RCW 89 near the end of the jet. 
Assuming that the synchrotron nebula surrounding the pulsar 
can be described in terms of a torus and jets similar to the Crab 
pulsar, the morphology suggests a large angle between pulsar 
spin axis and line of sight. 

The detection of pulsed emission at hard X-rays was first 
reported by Kawai et al. (1991) using Ginga LAC 2-60 keV 
data. They found that the X-ray pulse lags the radio pulse 
by 0.25 ± 0.02 in phase. The spectrum of the pulsed emis- 
sion could be represented by a power-law with photon index 



-1.3 ± 0.05 confirming its hard nature (ECawaietal. 1993) 



Wilson et al. (1993a,1993b) showed that pulsed emission was 
even detectable in the soft 7-ray regime using CGRO BATSE 



data (20-740 keV), confirmed later by Ulmer et al. (1993) and 



Matz etal. (1994) using CGRO OSSE data. The OSSE/BATSE 



pulse phase distribution (lightcurve) showed a phase offset 
with respect to the radio-pulse of 0.32 ± 0.02, slightly larger 
than the value obtained from the Ginga data. The OSSE spec- 
trum above 50 keV of the pulsed emission could be described 
by a power-law with photon index —1.68 ± 0.09, consis- 
tent with the spectral findings, a = —1.64 ± 0.42, from 
the balloon-borne Welcome instrument (94-240 keV) as re- 



ported by 3unji et al. (1994) The spectral measurements by 
EXOSAT, Ginga, Welcome, OSSE/BATSE suggest a spectral 
steepening (softening) towards higher energies. 

Recently, Rossi XTE absolute timing results on PSR 



B1509-58 were presented by Rots et al. (1998) showing high 



resolution pulse profiles for energies between 2 and 128 keV. 
A comparison of the profiles measured by RXTE (2-16 keV) 
and BATSE (> 32 keV) was made and a 0.03 phase shift 
of the BATSE lightcurve w.r.t. the RXTE lightcurve, peaking 



at 0.27 ± 0.01, was found by the authors. Rots etal. (1998) 



also performed a pulse-phase resolved spectral analysis using 
PCA and HEXTE data. The photon indices of the power-law 



fits to the 2 - 200 keV data in various 0.05 wide phase slices 
within the pulse are all consistent with one single value of 
a = —1.345 ± 0.01. This value is in line with the value of 



-1.3 ± 0.05 found by Kawai etal. (1991) RXTE spectral re- 



sults for both the pulsed and unpulsed component are described 



in detail by Marsden et al. (1998) the spectrum of the pulsed 



component (radio phase range 0.17 and 0.53) could be repre- 
sented by a power-law with index —1.358 ± 0.014 with no ev- 
idence for a spectral break seen up to ~ 200 keV, while the 
index of the unpulsed component was —2.215 ± 0.005. 

At high-energy gamma-rays (> 30 MeV) Brazier et al. 
1994 reported only upper-limits for pulsed emission using 
CGRO EGRET data from 3 viewing periods (VP 12, 23 and 
27; see e.g. Table |lj) during the all-sky survey of CGRO. A 
study by Fierro (1995) analysing EGRET data from Cycle I-III 
yielded an interesting ~ 4er source feature for energies above 
100 MeV consistent in position with the pulsar. Timing analy- 
ses resulted in non-detections for pulsed emission in the 30- 1 00 
MeV, > 100 MeV and > 1 GeV energy windows. 

The detection of hard pulsed emission below ~ 700 keV 
and the non-detection above 30 MeV makes PSR B1509-58 a 
very interesting candidate for COMPTEL, the Compton Tele- 
scope aboard CGRO and sensitive to photons with energies be- 
tween ~ 0.75 and 30 MeV. COMPTEL analysis using view- 
ing period 23 (see Table [j]) data only yielded an interesting 
J> 3cr detection of pulsed emission in the 0.75-1 MeV en- 
ergy window with a pulse aligned with the pulse observed by 
BATSE/OSSE ( formsen et al. 1994 |Carramiflana et al. 1995[ ). 
However, a timing analysis of COMPTEL 0.75-30 MeV data 
from observations spread over more than 4 years yielded only 
a marginal detection of pulsed emission at energies below 3 
MeV (Carraminana et al. 1997). Here, we will report on the 



COMPTEL timing- and spatial analyses using all available data 
up to and including CGRO Cycle-6 data. Prompted by our re- 
sults in the highest standard energy interval (10-30 MeV) in the 
COMPTEL analysis, we also (re)analysed all publicly available 
EGRET data on this source. 



2. Instrument description and observations 

COMPTEL is the imaging Compton Telescope aboard CGRO 
and is sensitive for 7-ray photons with energies in the 0.75- 
30 MeV range. Its detection principle relies on a two layer 
interaction: a Compton scatter in one of the 7 upper-detector 
(Dl) modules followed by a second interaction in one of the 
14 lower-detector (D2) modules. Main measured quantities are 
the angles (x, ip) specifying the direction of the scattered pho- 
ton (from the interaction loci in Dl and D2) and the energy de- 
posits in the D1/D2 modules where the interactions took place. 
From the last two quantities we can calculate the scatter angle 
Tp and the total energy deposit Etot (see for a full description 
|Schonfelder et al. 1993| ). 



Its energy resolution is 5-10% FWHM and due to its large 
field of view of typically 1 steradian it is possible to monitor 
a large part of the sky simultaneously with a position determi- 
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Table 1. COMPTEL observations used in current study with PSR B 1509-58 less than 30° off-axis 
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f TJD = JD - 2440000.5 = MJD - 40000 



* EGRET in narrow field mode; opening angle FoV is 19° 



nation accuracy of ~ 1°. The events are time tagged with an 
accuracy of 0. 125 ms. 

In this study we selected those CGRO Cycle I- VI view- 
ing periods for which the angle between the pointing axis (co- 
aligned with the COMPTEL/EGRET z-axis) and PSR B1509- 
58 is less than 30°. Details for each individual observation can 
be found in Table [j]. The last but one column specifies the ef- 
fective exposure in the 3-10 MeV energy window assuming a 
E~ 2 dependency of the photon flux. The calculation took into 
account Earth blocking effects and utilizes the timeline infor- 
mation. The last column indicates the status of the sparkcham- 
ber of the EGRET high-energy (30 MeV - 30 GeV) instrument 
aboard CGRO. EGRET data (with the sparkchamber ON) from 
the first 4 CGRO observation Cycli have been retrieved from 



the COMPTON Science Support Center and have subsequently 
been used in spatial- and timing analyses. 

3. COMPTEL timing analysis 

3.1. Event selections 

Prior to the actual timing analysis we have to specify the event 
selection criteria to which the events are subjected. The se- 
lection criteria applied here are the same as those in the tim- 
ing analysis of PSR B 195 1+32 (see Kuiper 1998a) except for 
one selection parameter, namely farm, the difference angle be- 
tween the calculated scatter angle Tp and the geometrical scat- 
ter angle ip geo . The last quantity can be determined from the 
known source position ( (J,b)=(320?321,-l?162)) and the scat- 
ter direction angles (x, ip). For a point-source the distribution 
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Table 2. PSR B 1509-58 radio-ephemerides 
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t Ephemeris has been provided by V. Kaspi (private communication). 

* Epheme ris has been been derived from radio timing data using the ATNF Parkes radio telescope, 
(see also http://www.atnf.csiro.au/researcli/pulsar/psr/archive/ ) 



The dispersion measure used in the calculation of absolute phases was 253.2 pc cm 3 . 



of (farm (i.e. the ARM-distribution) is a narrowly peaked distri- 
bution with a maximum near ip arm = and a wing for positive 
farm values due to incompletely absorbed events. The imag- 
ing capabilities of COMPTEL rely on this sharp asymmetric 
distribution of ip ar m- The relative contributions of the peak and 
wing, and the width of the peak are a function of input photon 
energy. This means that instead of fixing |v? ar m| to a value in 
the range 2? 5 to 3? 5 irrespective of the selected energies, as 
turned out to be the optimum range from COMPTEL studies 
on the Crab ( |Muchetal. 1995| ) and Vela ( |Kuiperet al. 1998h| ) 
pulsars, an energy window dependent ARM selection is more 
appropriate. In this study we have determined a priori the op- 
timal value of |<^ arm | for each energy window by estimating 
the maximum in the Signal-to-Noise vs. |</? arm | relation. The 
latter relation can be derived from a 3 dimensional (x>ip,¥>) 
point source model for the energy window involved and the 
total measured 3d-event distribution in the same energy win- 
dow, heavily dominated by instrumental background events 
(90-95%). The following energy dependent criteria on |</? orm | 
appeared to be appropriate: 3?5 for the energy window 0.75- 
1 MeV and 2?5 for the energy windows 1-3, 3-10 and 10-30 
MeV. The fraction of counts from a point-source within the 
ARM-cut is typically ~ 60%. It is also worth mentioning that 
the ARM cut applied in the timing analysis reduces the number 
of events handled in the timing analysis to typically 10% of the 
number of events available for the imaging or spatial analysis, 
in which the full 3d-dataspace is employed. 

For the 10-30 MeV interval we have departed from the 
"standard" Time of Flight (TOF) and Pulse Shape Discrim- 
ination (PS D) windows (see for a description of these event 
parameters |Schonfelder et al. 1993| ) of 113-130 and 0-110 re- 
spectively, normally applied in the timing analysis (see e.g. 
Kuiper 1998a) an d have used the optim um TOF and PSD win- 
dows derived by Collmaretal. (1997) in their study on opti- 
mum parameter cuts using the Crab pulsar/nebula signature in 
the COMPTEL event space. 



3.2. Pulse phase folding 

Once the event selection criteria were settled we proceeded as 
follows: the arrival times (recorded with an intrinsic resolution 
of 0.125 ms) of the events passing through our selection filters 



are converted to arrival times at the Solar System Barycentre 
(SSB) using the known instantaneous spacecraft position, the 
source position and the solar system ephemeris (JPL DE200 
Solar System Ephemeris). The pulse phase cf) is calculated from 
the following timing model: 



v ■ At + 1/2 ■ v ■ At 2 + 1/6 • v ■ At 3 



(1) 



In this formula At is given by At = t e — t° with t e the event 
SSB arrival time and t° the reference epoch. The values em- 
ployed here for t°, v, z>, i), <fio are given in Table |[ The RMS 
error of the timing models listed in Table || is typically 10 milli- 
periods or 1.5 ms, sufficiently accurate to keep coherency and 
allowing pulse phase folding over long time spans indicated by 
the validity range. 

3.3. Pulse profiles in the 0.75-30 MeV energy range 

The pulse phase distribution resulting from phase-folding 
COMPTEL Cycle I- VI 0.75-30 MeV data is shown in Fig. [jj 
The modulation significance of the unbinned sample of pulse- 
phases is 5.4c employing a Z^-test (Buccheri et al. 1983) with 



2 harmonics. This is the first firm detection of pulsed emission 
above 0.75 MeV from PSR B 1509-58. The puls e is roughly 



aligned with the pulse observed by OSSE/BATSE ( Ulmer et al 



1993) and peaks at phase 0.38±0.03 (obtained from a Gaussian 
+ background fit). 

We have split the integral energy window of 0.75-30 MeV 
up into 3 smaller energy windows, 0.75-3 MeV, 3-10 MeV 
and 10-30 MeV and performed similar timing analyses. The 
modulation significances (Zf -test) found for the 3 energy win- 
dows are 3.7<r, 4.0u and 2.1c respectively, proving the detec- 
tion of pulsed emission up to at least 10 MeV. The lightcurves 
are shown in Fig.|[ The 10-30 MeV lightcurve, having at face 
value a non-significant modulation, shows an indication for 
an enhancement in the phase range containing the pulse at 
lower energies. However, a narrower pulse might be visible 
near phase 0.85, which is absent at lower energies. A more de- 
tailed discussion on the 10-30 MeV lightcurve will be given 
once the spatial analysis, enabling the measurement of the to- 
tal flux from a source, has been introduced. 
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Fig. 1. Radio-aligned (radio pulse at phase 0) 0.75-30 MeV 
COMPTEL lightcurve of PSR B 1509-58. A double cycle is 
shown for clarity. Data from Cycle I- VI observations (see Table 
|lj) have been used. The modulation significance is 5.4<r adopt- 
ing a Z|-test. A typical error bar is indicated. 



3.4. Pulse profiles from soft X-rays to medium energy "f-rays 

Rots et al. (1998) studied the pulse shape of PSR B1509-58 
as a function of energy at medium and hard X-ray energies. 
He found that the RXTE 2-16 keV pulse peak lags the ra- 
dio peak by 0.27 ± 0.01 in phase. This value is consistent 
with the lag of 0.25 ± 0.02 foun d in Ginga 2-60 keV data by 
Kawai et al. (1991) . Ulmer et al. (1993) found at soft 7-ray en- 
ergies a phase lag of 0.32 ± 0.02 in CGRO BATSE and OSSE 
data. The difference of 0.05 ± 0.022 (2.3a effect) between 



these values was considered troublesome by Rots et al. (1998" . 
It could not be attributed to a CGRO clock absolute timing 
uncertainty and this triggered Rots et al. (1998) to reprocess 5 
years of BATSE data yielding now a slightly smaller phase lag 
of 0.30 consistent both with Ulmer's previous estimate and also 
with the RXTE estimate. 

The remaining (insignificant) difference of 0.03 ± 0.022 
(= 5 ms) still seems to be too large to be explained by 
CGRO/RXTE clock uncertainties and could be due to intrin- 



sic pulsar emission properties. Rots et al. (1998) also investi- 



gated the effect of the pulsar dispersion measure on the phase 
offset of the radio- and RXTE X-ray pulse. Using the most 
recent value of the dispersion measure of PSR B1509-58 of 
255.3 ± 0.3 pc cm~ 3 instead of the widely used value of 
253.2 ± 1.9 pc cm' 3 results in an offset of 0.29 ± 0.01, in- 
dicating that uncertainties in the dispersion measure can result 
in uncertainties in the radio/X-ray pulse phase lag of ~ 0.02. 

The RXTE 2-16 keV profile consists of a broad asymmet- 
ric pulse: a rapidly rising pulse reaching its maximum near 
phase 0.27 followed by a more gradual decline with a "shoul- 
der" near phase 0.40. This profile is shown in Fig.||b along with 
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Fig. 2. Radio-aligned COMPTEL Cycle I- VI lightcurves of 
PSR B 1509-58 for 3 different energy windows: 0.75-3 MeV 
(a), 3-10 MeV (b) and 10-30 MeV (c). Double cycli are shown 
for clarity, while typical error bars are indicated. The modula- 
tion significances adopting a Z|-test are 3.7er, 4. Oct and 2.1cr 
for the energy windows 0.75-3 MeV, 3-10 MeV and 10-30 
MeV respectively. This proves that pulsed emission is detected 
at least up to 10 MeV. 



the X-ray /soft 7-ray profiles measured by different instruments 
in different energy windows. It is clear from this collage that 
the ASCA 0.7-2.2 keV soft X-ray profile is in detail very sim- 
ilar to the RXTE 2-16 keV hard X-ray profile. The BATSE > 
32 keV soft 7-ray profile also resembles the RXTE profile ex- 
cept for an apparent overall offset of 0.03 with respect to the 
RXTE profile. The COMPTEL 0.75-10 MeV lightcurve on the 
other hand reaches its maximum near 0.38 which coincides in 
phase with the shoulder clearly visible in the RXTE lightcurve. 
Whether this apparent offset has a statistical origin (cf. the typ- 
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Table 3. PSR B 1509-58 pulsed fluxes derived from the timing 
analysis with the "pulsed" interval defined to be phases 0.15- 
0.65. 



_B-window 



Flux 

ph/cm 2 ■ s ■ MeV 



0.75 - 3 MeV 
3.0 - 10 MeV 
10.0 - 30 MeV 



(3.69 ±0.73) x 10~ 5 
(4.52 ± 0.77) x 10" 6 
(1.21 ± 0.85) x 10~ 7 



0.00 0.40 



0.80 1.20 1.60 2.00 



In order to study the difference in morphology between the 
COMPTEL 0.75-10 MeV and RXTE 2-16 keV lightcurves we 
have fitted the RXTE 2-16 keV profile in terms of a back- 
ground and two Gaussians (7 free parameters). This resulted 
in a narrow component peaking at phase 0.250 ± 0.008 with 
width 0.056 ±0.008 and a broader component at 0.386 ± 0.012 
and width 0.129 ± 0.006. The first narrow pulse accounts for 
25.7 ± 4.3% of the total pulsed emission. A similar fit has been 
performed on COMPTEL 0.75-10 MeV data, but now with po- 
sitions and widths fixed to the values found in the RXTE 2-16 
keV fit (3 free parameters). In this case the first narrow pulse 
can account for only 13 ± 18% of the total pulsed emission, 
consistent with being absent, and the profile can satisfactorily 
be described by just the broad second pulse near 0.39. This 
strongly suggests that the pulse shape changes from soft X-rays 
to medium energy 7-rays. 



Phase 

Fig. 3. Radio-aligned lightcurves of PSR B 1509-58 from soft 
X -rays to mediu m energy 7-rays: a) ASCA 0.7-2.2 keV 



dSaito et al. 1997] ; private communication), b) RXTE 2-16 keV 
(Rots et al. 1998), c) BATSE > 32 keV (Rots et al. 1998) and 
d) COMPTEL 0.75-10 MeV. The dashed vertical line indicates 
the phase lag of 0.27 found in the RXTE lightcurves. Note 
that the "center of mass" of the pulse appears to shift towards 
higher pulse phases for higher energies, with the maximum of 
the pulse at MeV-energies (COMPTEL) coinciding in phase 
with the "shoulder" visible in the pulse shape at X-ray ener- 
gies. Typical la error bars are indicated in each figure at phase 
0.95. 



ical error bar in the COMPTEL 0.75-10 MeV lightcurve which 
has a mean level of ~ 93000 counts in this 15 bin lightcurve) 
or is due to an intrinsic property of the pulsar's high-energy 
emission is difficult to decide. The intrinsic timing-resolution 
of 0. 125 ms of CGRO/COMPTEL is sufficiently accurate to al- 
low the lightcurve to be binned in several hundreds of bins and 
can not be responsible for the offset. Also, BATSE and COMP- 
TEL use the same CGRO clock. 



3.5. Pulsed 0.75-30 MeV fluxes from the excess counts in the 
lightcurves 

Based on the RXTE 2-16 keV lightcurve we have defined a 
"pulsed" and an "unpulsed" phase interval in the pulse phase 
distribution: the pulsed interval extends from phase 0.15 to 
0.65 and the unpulsed (background) interval from 0.65 to 1.15. 
This break-down is such that for a pulse shape as measured by 
RXTE 90.4% of the pulse is located in the "pulsed" interval. 

Applying this definition to the COMPTEL data we can de- 
termine the pulsed excess counts in various energy slices by es- 
timating the underlying background as the averaged level in the 
unpulsed part of the lightcurve. We derived these pulsed excess 
counts for the 0.75-3, 3-10 and 10-30 MeV energy windows 
and converted these to pulsed flux values taking into account 
efficiency correction factors due to the applied ARM cuts in 
the timing-analysis (see Table ||). 

The weak 10-30 MeV flux value should be treated with 
care because we do not detect significant modulation (2.1u). 
Moreover, the lightcurve shows indications for a second pulse 
in the "unpulsed"(background) phase interval. If this pulse is 
genuine, then the true flux is underestimated (see next section). 



4. COMPTEL spatial analysis 
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4.1. Analysis method 

The measured event parameters (x,ip,ip,Etot) constitute a 4- 
dimensional data space, in which we have to search for a 
"source signature". In practice the dimension of the data space 
is lowered by integrating along the Etot direction between 
user selected boundaries. The event distribution of a point 
source (the Point Spread Function, PSF) in this reduced 3-d 
data space (x,i>,Tp) is concentrated in a cone-shaped struc- 
ture with its apex at the source position (xo,ipo)- In the spa- 
tial or imaging analysis we proceed as follows: We gener- 
ate a background model from the sparsely filled event datas- 
pace (3d) through sophisticated smoothing techniques (see e.g. 
Bloemen et al. 1994 1. Because the measured events in the 3d- 
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dataspace are primarily internally generated background events 
(~ 90—95%) this background model represents a good approx- 
imation of the genuine instrumental background. The search 
for point sources in the measured 3d-dataspace is accomplished 
by a maximum likelihood ratio (MLR) test at scan positions in 
the selected sky field. In the null hypothesis Ho the measured 
events are described in terms of a background model alone, 
while in the alternative hypothesis Hi the data are described in 
terms of a background model and a point source at a given scan 
position. From the parameter optimizations under Hi we can 
derive the source flux. From the optimized likelihoods under 
Hi and Hq we can determine the maximum likelihood ratio 
A at each scan position, giving us information on the detec- 
tion significance of a source. For more detailed information see 
Kuiper etal. (1998a). 



4.1.1. Imaging results for energies < 10 MeV 

For consistency purposes we have produced MLR maps for 
the same energy windows (0.75-3 and 3-10 MeV) as used in 
the timing analysis, and also made selections on pulse phase, 
the "pulsed" and "unpulsed" windows introduced in Sect. 3.5, 
while we also considered the total emission, the sum of both. 
In Fig. |] the MLR images are shown for the "pulsed" and "un- 
pulsed" phase intervals in the energy windows 0.75-3 and 3-10 
MeV. The contours start at an equivalent of 3a in steps of la 
for 1 degree of freedom (i.e. at A — 9, 16, 25, 36 etc.), repre- 
sentative if the source position is a priori known. 

From these maps it is evident that significant emission at 
the pulsar position (indicated by a " X " mark) is only observed 
in the "pulsed" maps: 4.5<r in the 0.75-3 MeV energy window 
and 8a in the 3-10 MeV window. From the maximum likeli- 
hood fits we can also obtain estimates for the source flux. How- 
ever, in the maps shown in Fig. ^ the structured galactic diffuse 
emission is not included in the background model, and may 
contribute to the source flux. In addition, any pulsar/nebula 
DC-emission will contribute equally to the "pulsed" and "un- 
pulsed" maps. Since A is not zero at the source position in the 
"unpulsed" maps, we subtracted the measured fluxes (counts) 
in the "unpulsed" maps from those in the "pulsed" maps in or- 
der to obtain independent estimates of the "pulsed" fluxes. For 
the 0.75-3 MeV interval we found 7202 ± 1386 counts corre- 
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Fig. 4. COMPTEL "pulsed" (phases 0.15 - 0.65) and "un- 
pulsed" (0.65 - 1.15) MLR maps for energies below 10 MeV. 
The upper panel shows the results for the 0.75-3 MeV energy 
window, the lower panel the 3-10 MeV analogon. The left side 
shows the pulsed maps and the right side the unpulsed maps. 
The solid contours start at the 3cr-significance level (A = 9, 1 
degree of freedom) in steps of la. Significant emission at the 
pulsar position (x symbol) is only observed in the "pulsed" 
maps. 



lating with a point source at the pulsar position in the "pulsed" 
interval and 2264 ± 1382 in the "unpulsed" interval. The dif- 
ference of 4938 ± 1960 is consistent with the pulsed excess 
counts in the timing analysis. The number of counts found in 
the "unpulsed" window can be fully attributed to contributions 
of galactic diffuse origin, as has been verified in a simulation 
of the galactic diffuse emi ssion using the model components 
and scale factors found by Bloemen et al. (1999) This means 
that there is no room for any significant pulsar or nebula DC- 
component in this energy window. 

For the 3-10 MeV energy interval we obtained a similar 
picture. In this case we find 5233 ± 593 counts correlating with 
a point source at the pulsar position in the "pulsed" map and 
2755 ± 590 in the "unpulsed" map. This leaves 2478 ± 840 
counts for the genuinely pulsed emission, consistent with the 
timing analysis results. The "unpulsed" value can again be ex- 
plained with being of galactic diffuse origin. 

In conclusion, the fluxes derived from the spatial analysis 
for energies below 10 MeV are consistent with those derived 
from the timing analysis and there is no evidence for the detec- 
tion of DC-emission from the pulsar or its nebula. 
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Fig. 5. COMPTEL "total" MLR map for the energy win- 
dow 10-30 MeV, presenting the source detection significance 
above the instrumental and galactic diffuse background. A 5.6er 
source is detected at the pulsar position. The solid contours start 
at the 3r7-significance level (A = 9, 1 degree of freedom) in 
steps of la. The pulsar position is marked by the x symbol. 



4.1.2. Imaging results for energies > 10 MeV 

In the 10-30 MeV energy window we did not detect a signifi- 
cant (2.1a) modulation in the light curve, but an enhancement 
is visible in the pulse phase interval in which the pulse is con- 
centrated at lower energies, and in the "unpulsed" interval a 
high bin shows up near phase 0.85. A MLR image (total) in this 
energy window yielded somewhat surprisingly a strong source 
feature (~ 6er) consistent in position with the pulsar, by sim- 
ply fitting the source contribution above the instrumental back- 
ground model. 

It should be noted here that the 10-30 MeV energy win- 
dow is COMPTEL's "cleanest" window not polluted by time- 
varying instrumental background lines contrary to the energy 
intervals below 10 MeV. Estimates for the galactic diffuse 
emission in this en ergy window can be considered reliable (e.g. 
Strong et al. 1994 , Bloemen et al. 1999) and are consistent with 
extrapolations towards lower energies of EGRET > 30 MeV 
measurements. 

When we now include also the total galactic diffuse back- 
ground model in the total background, the source remains (see 
Fig. |]). At the pulsar position a detection significance of 5.6<r 
is reached. A break-down into "pulsed" and "unpulsed" maps 
yields sources at the pulsar position with statistically equal 
fluxes in both maps. This can mean that we either detected DC- 
emission from the pulsar and/or its surrounding synchrotron 
nebula or emission from a nearby unrelated source. But, in 
these cases we would expect to see emission below 10 MeV 



where COMPTEL is more sensitive. Another possibility is that 
we actually detected pulsed emission but with a different pulse 
profile. 

We investigated the latter possibility further. In this respect 
it is instructive to go back to the 10-30 MeV lightcurve and 
superpose the background level as determined in the spatial 
analysis from the total number of counts correlating with a 
source at the pulsar position above the instrumental and galac- 
tic diffuse background. The lightcurve with this spatially deter- 
mined background level (and its ler-errors) is shown in Fig. |8^. 
The existence of the main pulse, definitely detected below 10 
MeV, is again very suggestive for energies above 10 MeV. The 
number of excess counts in the "pulsed" interval (0.15-0.65) 
above of this spatially determined background level is 279 ± 57 
counts, resulting in a flux of (3.37 ± 0.70) • 10 -7 ph/cm 2 ■ 
s ■ MeV. In the "unpulsed" region we find then a flux of 
(2.16 ± 0.70) • 10~ 7 ph/cm 2 ■ s ■ MeV, predominantly due 
to the excess near phase 0.85. This ~ 3er feature above the spa- 
tially determined background, only shows up above 10 MeV 
and, if genuine, would indicate a new pulse component with a 
completely different spectral behaviour from that of the main 
pulse. This peculiar behaviour of the pulsed emission above 10 
MeV in the COMPTEL energy range motivated us to analyze 
also the contemporaneous > 30 MeV EGRET data. 



5. Analysis of EGRET data 

No significant (pulsed) emission above 30 MeV has been re- 



ported from EGRET data, but, two papers ( Brazier et al. 1994 1 
^ierro 1995 ) report on weak (irregular) source features (~ 3<t) 



near PSR B 1509-58 for energies above 100 MeV. It is worth 
mentioning that in the first paper the analysis was constrained 
to data from CGRO Cycle I, namely VP's 12, 23 and 27 (cf. 
Table 1), while in the second paper data from Cycles I to III 
were analyzed. Here, we reanalyzed all available Cycle I-IV 
EGRET data, the maximal exposure on the source. The obser- 
vations used are those in Table 1 in which the spark chamber 
high-voltage was enabled (see last column Table 1 ; the data of 
Cycle V do not contribute since EGRET operated in its narrow 
field mode with an effective FoV opening angle of 19°). Stan- 
dard event selections were applied, however, the requirement of 
a minimum energy deposit of 6.5 MeV in the TASC was aban- 
doned (for a description of EGRET see Thompson et al. 1993 ). 
This selection criterion is mainly effective for background sup- 
pression for measured total energies above 100 MeV. We ver- 
ified using the Crab pulsar that for energies below 100 MeV 
the detection of the pulsed signal significantly increases when 
the TASC threshold of 6.5 MeV is ignored. In order to suppress 
Earth albedo 7-rays, standard Earth horizon angle cuts were ap- 
plied for the differential energy ranges 30-50 MeV, 50-70 MeV 
and 70-100 MeV, which are roughly equivalent to a 3<r cut. For 
consistency, the same selections were applied in the spatial and 
timing analyses. 
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5.1. Spatial analysis 

The imaging analysis of EGRET data is similar to that of 
COMPTEL data with the main difference that the EGRET anal- 
ysis is performed in a 2-d dataspace. After reformatting the 
EGRET data and response for import i n the COMPTE L Anal- 
ysis and Software System COMPASS ( |de Vries 1994] ) we can 
use the same analysis programmes. One important difference 
is that the EGRET data are almost free of instrumental back- 
ground. This means that the point source emission is searched 
for above of the dominating galactic and extragalactic 7-ray 
backgrounds only. The spatial structure of the first component, 
predominantly due to the interaction of cosmic rays and inter- 
stellar Hydrogen, is approximated by the measured total col- 
umn density of atomic Hydrogen Hi and molecular Hydro- 
gen as traced by the CO-molecule, identical to our analysis 
of COMPTEL data. The latter component can simply be de- 
scribed by an isotropic model. 

The MLR map for the 30-100 MeV energy window, ad- 
jacent to the COMPTEL 10-30 MeV window, combining all 
available Cycle I-IV data is presented in Fig. ||. A strong 6.7a 
excess shows up near PSR B1509-58. Apart from the pulsar a 
few other identified and unidentified EGRET sources seem to 
be visible in the map (see figure caption). 

The most likely sources responsible for the excess near 
(I, b) = (320, 0) are PSR B1509-58 and 2EG J1443-6040 at 
(I, b) = (316.28, —0.75). Fitting both sources simultaneously 
above of the galactic and extragalactic backgrounds, together 
with emissions from Cen-A and 2EGS J 1429-4224 yields a 
4.4cr excess at the position of PSR B1509-58 with a flux of 
(3.3 ± 1.0) • 10~ 8 ph/cm 2 ■ s ■ MeV. This is shown in the 
MLR map of Fig. [7] in which the other three fitted sources and 
background emissions are "subtracted". The detection signifi- 
cance of the source correlating with 2EG J1443-6040 is 2.6er 
and has a flux of (1.9 ± 1.0) • 10~ 8 ph/cm 2 ■ s ■ MeV. It 
should be noted that this unidentified EGRET source was only 
detectable in VP's 12 and 27 and not during observations per- 
formed later on, indicating its likely transient behaviour. In fact, 
it does not appear anymore in the third EGRET source catalog 
( Hartman et al. 1999). (Note: After submission of this paper, 
the EGRET group, using analysis based on the third EGRET 
catalog, confirmed the finding of a 4.3(7 significance source in 
the 30-100 MeV band, consistent in position with PSR B 1509- 
58; D.J. Thompson, private communication). 

We performed a similar imaging study in the 100-300 MeV 
energy window and found also a ~ 2a source feature near PSR 
B 1509-58 above of the galactic and extra-galactic background, 
as well as the previo usly reported EGRE T s ources, consi stent 
with the findings of Brazier et al. 1994 and Fierro 1995 . The 
100-300 MeV flux estimate for PSR B1509-58 is (5.9 ± 3.5) • 
10~ 10 ph/cm 2 ■ s ■ MeV. Even in the 300-1000 MeV window 
a marginal excess is visible consistent in position with PSR 
B1509-58 (flux (4.5 ± 3.6) • 10~ n ph/cm 2 ■ s ■ MeV). Above 
1000 MeV no excess is found. These low flux values can also 
be converted in upper-limits. 




Fig. 6. EGRET "total" MLR map for the energy window 30- 
100 MeV, showing the likelihood of a source detection above 
of the galactic and extragalactic diffuse backgrounds. Contours 
of A = 9, 16, 25 (= 3, 4, 5 a); Solid contours positive values, 
broken ones negative values. Several identified and unidenti- 
fied sources are indicated by different symbols: PSR B1509- 
58 (x), Cen-A (A), and the unidentified EGRET sources 
2EG J1412-621 1 (□), 2EG J1443-6040 (+), 2EGS J1418-6049 



(o) and 2EGS J 1429-4 224 (*) (see [Thompson et al. 1995| and 
[Thompson et al. 1996a ). A 6.7<r source feature near (l,b) = 
(320, 0) shows up and is likely composed of the contributions 
of PSR B1509-58 and the unidentified EGRET source 2EG 
J 1443-6040. 



5.2. Timing analysis 

Given the detection of a significant source in the spatial anal- 
ysis for energies between 30 and 100 MeV, we first selected 
events in this energy range to search for the modulation. The 
event selection procedures in the EGRET timing analysis are 
equivalent to those applied in the COMPTEL timing anal- 
ysis except that we are now dealing with a 2-dimensional 
dataspace. We applied an energy dependent cone selection 



( Thompson et al. 19961: ) roughly selecting 68% of the source 



counts. In this 30-100 MeV energy interval we do not find 
significant modulation (l.lcr), although the shape of the pulse 
phase distribution bears some resemblance to the COMPTEL 
10-30 MeV lightcurve. Both are shown in Fig. ^. Estimates 
of the underlying backgrounds as determined from the spa- 
tial analysis fitting PSR B1509-58 and 2EG J1443-6040 si- 
multaneously are indicated. Also shown in the figure is the 
combined 10-100 MeV lightcurve, which has a modulation 
significance of 2.3a for 2 harmonics in the Z^-test. Again, 
the summed background level determined in the spatial anal- 
yses is indicated, as well as the Kernel Density Estimator 
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Fig. 7. EGRET MLR map for the energy window 30-100 MeV, 
as in Fig. g. The contributions from 2EG J 1443-6040, Cen-A 
and 2EGS J 1429-4224 are subtracted. At the position of PSR 
B1509-58 (x) a significance for a source detection of 4.4<r re- 
mains. 



(KDE, de Jager et al. 1986) of the unbinned pulse-phase dis- 
tribution with the ±lcr error bounds. The KDE provides an 
(asymptotically) unbiased view on the genuine pulse shape. 
This lightcurve suggests a double-peaked profile: an enhance- 
ment in the previously defined "pulsed" interval, and a pulse 
near phase 0.85, which has become more pronounced in this 
summed lightcurve. 

Even though the modulation significance of this 10-100 
MeV lightcurve is still marginally significant, the apparent 
double-peaked structure above the background level estimated 
in the spatial analysis, makes it very suggestive that the source 
detected in the spatial analysis is PSR B1509-58 with pulsed 
emission upto the EGRET energies. At least, the spatial and 
timing analyses are consistent with this interpretation. If cor- 
rect, the lightcurve morphology changed from one broad single 
pulse to a profile with an additional pulse near phase 0.85. 

6. Spectral analysis 

Energy spectra can be derived from the spatial analysis as well 
as from the timing analysis. Given that there is some ambiguity 
in the interpretation of the results on PSR B 1509-58 for ener- 
gies above 10 MeV, we consider three cases: 

/) The COMPTEL "pulsed" spectrum of PSR B1509-58 
(0.75-30 MeV) determined by the excess counts in the phase 
window 0.15-0.65 above the average level in the complemen- 
tary phase interval. This represents the spectrum of the broad 
single pulse, which was already determined in section 3.5 (Ta- 
ble 3). In Fig. this spectrum is shown together with the spectra 
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Phase 

Fig. 8. a COMPTEL 10-30 MeV Cycle I- VI lightcurve (2.1cr): 
the background level and its ±lcr error estimated from the spa- 
tial analysis are indicated as a solid straight line and dashed 
lines, b EGRET 30-100 MeV Cycle I-IV lightcurve (1.1a): 
background level (and ±lcr errors) are determined in a spa- 
tial analysis including point sources at the pulsar and 2EG 
J1443-6040 positions, c Combined COMPTEL 10-30MeVand 
EGRET 30-100 MeV lightcurve (2.3a or 1.9% chance proba- 
bility). The summed background level from the spatial analyses 
is again indicated together with the smoothed curves represent- 
ing the KDE and its ±lcr error region. Typical error bars are 
indicated. 



(power-law fits) derived by other instruments from 1 keV upto 
~ 700 keV. Above 30 MeV are included the EGRET 2a upper- 
limits (D.Thompson-private communication) derived from spa- 
tial analyses using a subset of the observations listed in Table 
1 . The new COMPTEL flux values are consistent with an ex- 
trapolation of the OSSE power-law fit (PL-index of -1.68) upto 
10 MeV, followed by a drastic break, which will be discussed 
in section 7. 

ii) If the source detected by COMPTEL between 10 and 
30 MeV is PSR B 1509-58, then we have underestimated the 
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ROO collaboration in their spatial analysis for energies above 
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Fig. 9. The pulsed high-energy spectrum of PSR B1509-58 
from soft X-rays to hard 7-rays. The filled squares are the 
COMPTEL flux points as derived from the excess counts in 
the 0.15 — 0.65 phase range (case z), from section 3.5), while 
the open square represents the 10-30 MeV flux in the 0.15- 
0.65 phase interval above the spatially determined background 
(case if). The filled triangles are the 2a upper limits for the to- 
tal fluxes in the EGRET energy domain (D. Thompson - private 
communication). The polygons represents the best fit la error 
regions as measured by different instruments below the COMP - 
TEL energy win dow: ASCA (0.7- 10 keV; |Saito et al. 1997| ); 
Ginga (2-60 keV; fcawai et al. 1993b ; OSSE (50-750 keV; Matz 
et al. 1994); WELCOME (94-240 keV; Gunji et al. 1994); 
RXTE (2-250 keV; Marsden et al. 1998). The softening of the 
spectrum from soft X-rays to MeV gamma rays is evident, as 
well as the spectral break above 10 MeV. 



flux of the single broad pulse in the 10-30 MeV interval. De- 
termining then the flux (phases 0.15-0.65) above the spatially 
determined background (see Fig. ^) between 10 and 30 MeV, 
and inserting that in Fig. ^ gives a spectrum for the single broad 
pulse with a break energy between 10 and 30 MeV. In this case 
this spectrum does not represent the total pulsed spectrum for 
energies above 10 MeV. 

Hi) If the sources detected at the position of the pulsar by 
COMPTEL and EGRET in the spatial analyses above 10 MeV 
and 30 MeV, respectively, can be identified with the pulsar, then 
the total spectrum of PSR B1509-58 above 0.75 MeV can be 
derived by combining the COMPTEL flux values from the tim- 
ing analysis below 10 MeV (consistent with 100% pulsed) with 
the flux values from the spatial analysis of COMPTEL 10-30 
MeV data (section 4. 1 .2) and EGRET data above 30 MeV (sec- 
tion 5.1). This is shown in Fig. |l0| together with the "pulsed" 
spectra measured at lower energies by other instruments as well 
as the flux estimate from the detection (4. Icr) by the CANGA- 



1.5 TeV ( |Sako 1998| ). The increased flux level above 10 MeV 
in comparison with the spectrum shown in Fig. ^|is then mainly 
due to an additional pulse component near phase 0.85, having 
a completely different spectrum than the main pulse. For com- 
parison is also shown the "unpulsed" spectrum of PSR B1509- 
58 measured by RXTE upto 240 keV (Marsden et al. 1998). It 
is evident that, if the "unpulsed" emission extrapolates accord- 
ing to the RXTE spectrum towards the 7-ray regime, the 10- 
30 MeV flux measured at the PSR B1509-58 position can not 
be the (nebula) DC emission, supporting the pulsar interpreta- 
tion. Furthermore, the EGRET flux values above 100 MeV are 
consistent with the extrapolation of the fit to this "unpulsed" 
spectrum, which is even consistent with an extrapolation up to 
the claimed detection of the nebula by the CANGAROO col- 
laboration at TeV energies. Therefore, the spectral compilation 
in Fig. [l(] suggests that the combined COMPTEL / EGRET 
spectrum of PSR B 1509-58 represents the transition from dom- 
inantly pulsed emission below ~ 30 MeV to dominantly DC- 
emission above ^100 MeV. 



7. Summary and discussion 

The major findings based on primarily COMPTEL and to a 
lesser extent EGRET 7-ray data presented in this paper can be 
summarized as follows: 

I Pulsed 7-ray emission from PSR B 1509-58 has been de- 
tected with high significance up to 10 MeV as a broad 
asymmetric single pulse located in the phase interval where 
also the pulse occurs for X-ray and soft 7-ray energies up 
to ~ 700 keV. The pulse measured by COMPTEL between 
0.75 and 10 MeV reaches its maximum near radio phase 
0.38, shifted w.r.t. the value 0.30 measured by Rots et al. 
(1998) in the BATSE data above 32 keV. The COMPTEL 
maximum appears to coincide in phase wih the "shoulder" 
clearly visible in the RXTE 2-16 keV data. The pulsed 
spectrum is consistent with a power-law extrapolation up 
to lOMeV of the OSSE spectrum (PL-index of -1.68) mea- 
sured at lower 7-ray energies. 

II In the 10-30 MeV energy range we do detect a significant 
(5.6er) source positionally consistent with PSR B1509-58, 
however, the timing analysis yields a modulation signifi- 
cance of ~ 2a only. Based on just the COMPTEL data we 
cannot discriminate with certainty between the following 
interpretations: /) DC emission from the pulsar or its neb- 
ula, ii) emission from a nearby source, or Hi) pulsed emis- 
sion: the light curve exhibits an excess in phase with the 
pulse measured below 10 MeV and there is an indication 
for a second narrow component near phase 0.85. 

Ill In the adjacent 30-100 MeV EGRET energy band we 
detect a 6.7a source positionally consistent with PSR 
B 1509-58 which can be explained with contributions from 
PSR B1509-58 (4.4cr) and the nearby unidentified EGRET 
source 2EG J1443-6040 (2.6a). Timing analysis of the 30- 
100 MeV events yields a modulation significance of 1.1a 
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Fig. 10. The total COMPTEL and EGRET high-energy spectrum of PSR B 1509-58 (> 0.75 MeV, case in) including the spectral 
measurements below 0.75 MeV of the pulsed emission by other high-energy experiments (see caption Fig. ^|). Furthermore, the 
"unpulsed" spectrum of PSR B1509-58 measured by RXTE (Marsden et al. 1998) is shown, as well as the > 1.5 TeV DC-flux 
measurement by the CANGAROO collaboration (assumed to have a photon index of -2.5). The total emission breaks now around 
30 MeV due to the contribution from an additional pulse near phase 0.85 for energies in excess of 10 MeV. 



only. The combined COMPTEL/EGRET 10-100 MeV 
lightcurve (modulation significance 2.3<r) shows also the 
double-peaked shape: again a main pulse coinciding with 
the pulse observed at lower energies and a weak narrower 
pulse near phase 0.85. The results of the spatial and timing 
analyses are consistent with PSR B1509-58 being detected 
between 10 and 100 MeV, with a new pulse which is only 
visible between 10 and 100 MeV. 
IV The pulsed spectrum of PSR B 1509-58 shows a sharp break 
between 10 and 30 MeV. The break energy is close to 10 
MeV for the broad main pulse, and shifts to ~ 30 MeV 
for the total spectrum if the second narrow component is 
genuine. 

The implications of these new findings at medium and high- 
energy 7-ray energies will now be discussed. Theoretical mod- 
els explaining the high-energy electro-magnetic radiation from 
highly magnetized rotating neutron stars can be distinguished 
in two different catagories: 

a Polar Cap models (PC) in which the acceleration of charged 
particles along the open magnetic field lines in the vicinity 
of the magnetic pole(s) and subsequent cascades through 
high-energy radiation processes gives rise to the emerging 
7-ray spectrum. 

b Outer Gap models (OG) in which the acceleration of 
charged particles and the production of 7-rays takes place 
in vacuum gaps between the null-charge surface defined by 
VI ■ B = 0, with B the local magnetic field and O the pul- 



sar spin vector, and the light cylinder (R[ c = c/fi) along 
the last closed field lines. 

Both models rely on charge replenishment of the magneto- 
sphere mainly through e pairs. The polar cap models can be 
subdivided further primarily based on the energy reached by 
the particles (primaries) in the acceleration process. For pri- 
mary particle energies V <; 10 6 , with V the Lorentz factor 
or equivalently the dimensionless energy, energy loss in the 
form of high-energy 7-rays through curvature radiation is the 
most important radiation mechanism (CRPC models). For en- 
ergies 10 4 T < 10 6 energy loss through inverse Compton 
scattering with either thermal X-ray photons from the polar 
cap or from non-thermal cascade processes in the magneto- 
sphere will be more important (ICPC-models). Finally, for en- 
ergies 50 ^ r 10 resonant Compton scattering will be 
the dominant energy loss mechanism. Irrespective of the un- 
derlying energy loss mechanism of the primary particles it is 
assumed that the emergent 7-ray spectrum is softened by mag- 
netic pair production (7 —> e + + e~; Sturrock process) in the 
strong magnetic fields present in the vicinity of the magnetic 
poles. These pairs (secondaries) are produced in excited Lan- 
dau states, which decay by the emission of synchrotron photons 
which in turn can produce e ± pairs and a cascade can develop 
softening the initial 7-ray spectrum at each generation. 

Daugherty & Harding (1996) assume in a recent version 
of their CRPC model that a distant observer sees the emis- 
sion from just one polar cap. They also assume that the ac- 
celeration of the primaries occurs over an extended distance 
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above the polar cap surface and that the magnetic inclina- 
tion a (angle between ft and magnetic moment /x) is com- 
parable to the 7-ray beam opening angle 8 1 . The last two as- 
sumptions relax the requirement of both a m 8 1 and a small 
a, i.e. a nearly aligned rotator, in earlier versions of their 



For PSR B 1509-58 this CRPC model can not explain the 
observed low spectral cutoff energy occuring between 10 and 
30 MeV. However, the (very) strong polar surface magnetic 
field strength of ~ 3.1 x 10 13 GauB likely triggers another more 
exotic attenuation process to be active in the vicinity of the 
model ( paugherty & Harding 1994| ) to explain the 7-ray emis- ma g ne tic pole, namely photon splitting 7 -^7/ + 7", besides 



sion properties of pulsars 

Because the curvature radii R cr of the open field lines, 
assuming a dipole configuration, originating from the polar 
cap rim near 8 pc ~ axcsin(\/{R ns / Ri c )) are smaller than 
those originating at smaller polar angles (at the magnetic pole 
R cr — > 00 and no curvature emission can be produced) the 
initial CR 7-ray emission spectrum (oc 1/R cr ) is most pro- 
nounced and hardest (maximum of CR-spectrum is reached at 
w max — 0.29 • uj c oc r 3 /i? cr ) near the polar rim resulting in a 
hollow cone emission pattern. The mostly developed cascades 
originate near the polar rim which in turn soften the input CR- 
spectrum most efficiently. 

This means that a distant observer not only sees a double 
peaked profile when the hollow 7-ray emission cone passes his 
line of sight, but also a soft-hard-soft variation of the emission 
is expected, because the CR-emission emanating nearer the 
magnetic pole is less softened by cascading. Grazing through 
the 7-ray emission cone will result in a single pulse profile 
highly softened by cascade processes. The last situation may be 
applicable for PSR B 1509-58 requiring clear constraints on the 
geometrical aspects of both the system and the observer. In par- 
ticular, it demands that the pulsar spin axis / line-of-sight angle 
£ should approximately be equal to the sum of the magnetic 
inclination angle a and the 7-ray cone semi-angle 1 (beam 
radius). However, the large duty cycle of the single pulse of 
^ 0.4 poses an additional constraint, namely £ should then be 
smaller than a + 8 1 . This can be estimated using the following 
equation: 

8 1 = arccos(cos j3 — 2 sin a sin(a + /3) sin 2 (^-)) (2) 

which can be derived in the rotating vector model (e.g. Rankin 
1993). A$ is the duty cycle of the pulse, while j3 is the impact 
parameter defined as j3 = £ — a. Radio-polarization data yield 
in the context of this rotating vector model (Crawford et al. 
1997; private communication 1999) a ~ 18^3 degree and 
P ~ 20^20 degree. If we use the best fit parameters, then 7 ~ 
36°. Eventhough the uncertainties are large, the emitting rim 
must be very broad to explain the broad single pulse. 

The CRPC model also predicts a spectral cutoff in the 7-ray 
spectrum at several GeV. Another interesting feature is that a 
lower cutoff energy is expected for those pulsars with a higher 
magnetic field, because the softening by cascade processes is 
more efficient in these cases. Some of these general trends are 



attenuation by magnetic pair production. Harding et al. (1997) 



indeed observed in some 7-ray pulsars (Thompson et al. 1997) 



The current implementation of geometrical and physical as- 
pects, however, predicts lightcurves with a high degree of sym- 
metry. This symmetr y is a problem for e. g. the Vela pulsar high - 
energy lightcurve ( Grenier et al. 1988 ; fcanbach et al. 1994 ) 
which shows rather asymmetric bridge emission. 



showed that the photon splitting attenuation lengths can be 
shorter than the attenuation lengths for pair production for 
magnetic field strenghts B ^ 0.3B cr with B cr = m e c 3 /eh = 
4.413 x 10 13 G. This means that photon splitting acts as the 
dominant attenuation process for such strong magnetic fields 
and can suppress the emission of high-energy photons. Con- 
trary to the attenuation by magnetic pair production, photon 
splitting has no energy threshold and can degrade the pho- 
ton energy also below the magnetic pair production thresh- 
old (ut > 2m e c 2 / s'm(8kB), ®kB is the angle between pho- 
ton propagation direction and local magnetic field). Depend- 
ing on the splitting mode, partial and full splitting cascades 
are addressed, Harding et al. (1997)| calculate the high-energy 
spectrum of PSR B 1509-58 for various values of the model 
parameters 8kB and 8, the magnetic co-latitude angle, assum- 
ing that the initial high-energy photon emission originates from 
the neutron star surface in a polar rim at co-latitude 8. Our new 
medium-energy 7-ray data and in particular the spectrum of the 
main pulse severely constrain the magnetic co-latitude of the 
emission rim, irrespective of the splitting mode. A co-latitude 
of ~ 2°, close to the "classical" radius of the polar cap 8 pc , ap- 
pears to be required in the model calculations to be consistent 
with the evident spectral break between 10 and 30 MeV in the 
combined COMPTEL-EGRET spectrum. 

No detailed model calculations have been performed for 
PSR B1509-5 8 in the inverse Compton induced Polar Cap cas- 
cade scenario (Sturner & Dermer 1994; Sturner et al. 1995 and 
[Sturner 1995 ), although some interesting qualitative statements 
have been made. In particular, if B > 10 13 G and the neutron 
star surface temperature T > 3 • 10 6 K the Lorentz factors of 
the electrons are limited to 10 3 , which might explain the low 
cutoff energy in the spectrum of PSR B 1509-58 ( {Sturner 1995| ). 
However, this type of PC-model also suffers from predicting 
too symmetric lightcurves. The model requirement of a nearly 
aligned rotator (a ~ 8 1 ^ 5°) can formally not be excluded 
for the estimates of the magnetic inclination a ~ 18^]J degree 



(Crawford et al. 1997) and 0. 



36^30 degree. 



In outer gap scenarios (see e.g. Cheng et al. 1986a, 
pieng et al. 1986b| and |Ho 1989[ ) it is believed that stable vac- 
uum gaps (Holloway 1973) can be formed in the outer mag- 



netosphere along the boundary of the last closed field lines 
which extend from the null-charge surface to the light cylin- 
der. Voltage drops of typically ~ 10 13 V can be obtained 
across the gaps and accelerate e , created either in the gap 
or flowing in from beyond the light cylinder and from across 
the null-charge surface, to energies limited by curvature radia- 
tion and to a lesser extent inverse Compton scattering (off the 
ambient bath of lower energy photons). A geometrical calcula- 
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tion of the high-energy emission (beamed along the local mag- 
netic field in the outer gap resulting in a "fan" beam) from an 
outer gap region was successfull in reproducing qualitatively 
the observed 7-ray lightcurves of the known 7-ray pulsars 
( |Chiang & Romani 1994 ). In a recent paper Romani (1996) 



modelled the emergent high-energy emission based on curva- 
ture radiation reaction limited accelerated charges in the outer 
magnetosphere. In his model the primary particles (e^) emit 
curvature radiation with spectral cutoffs in the 1-10 GeV range. 
In order to tap the potential drop a small fraction (proportional 
to the optical depth t 77 and of the order of ^ 10~ 3 -10 -2 ) 
of these curvature photons interact with low-energy photons, 
provided by thermal surface emission from the neutron star or 
produced in the gap itself by synchrotron emission processes, 
to produce e ± pairs. The perpendicular momentum of the pro- 
duced pairs with respect to the local magnetic field is emitted 
in the form of synchrotron radiation peaking in the 1-10 MeV 
range. Part of the produced synchrotron radiation can inverse 
Compton scatter (ICS) off the primary particles to produce a 
TeV pulsed emission component (typically SS 1% of the pulsed 
GeV flux). The composite of the various emission components 
constitute the high-energy pulsed spectrum. Although no quan- 
titative estimates are presented, Romani (1996) argued that for 
short-period high-magnetic field pulsars like PSR B 1509-58 
the synchrotron flux produced in the gap itself will dominate 
the thermal surface emission increasing t 77 towards ~ 1 and 
thus significantly suppressing the GeV curvature component. 
This results in a synchrotron type spectrum peaking at MeV 
energies. From a geometrical point of view, the derived large 
values for the magnetic inclination and line-of-sight / pulsar 
spin axis angles seem to be naturally explained with an outer 
magnetospheric origin for the high-energy emission. 

The above considerations for both the polar cap and outer 
gap scenarios show that further model calculations are required 
to make decisive statements on the production site and mech- 
anisms of the high-energy pulsed radiation. In the small sam- 
ple of detected high-energy pulsars PSR B 1509-58 represents 
clearly a special case. This might now even be more so, given 
that in the above discussions we did not yet introduce the pos- 
sible detection of an additional pulse in the lightcurve for ener- 
gies above 10 MeV. This component must have a significantly 
harder spectrum than that of the broad main pulse. Such an 
additional narrow component may at first sight be more eas- 
ily incorporated in an outer-gap scenario. In this respect it may 
be of interest to refer to earlier claims of detection of multiple 
and variable components in the lightcurve of PSR B1509-58 
at TeV energies (> 1.5 TeV), which have been interpreted in 



the frame work of the outer gap scenario (de Jager et al. 1988, 



Nel et al. 1990). However, the earlier mentioned TeV observa- 



tions by the CANGAROO collaboration do not confirm these 
claims, and provide a significantly lower upper-limit to the 
pulsed emission. 

A long exposure of the source in the transition region of 
COMPTEL and EGRET is required to better study the change 
in morphology of the lightcurve. Additional COMPTEL obser- 
vations of this source have been scheduled (CGRO Cycle 8) 



and are aimed to provide better insights in the source charac- 
teristics especially above 10 MeV. 
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